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We study the heavy baryon mass spectrum on gauge configurations that include 2+1 flavors of 
dynamical improved staggered quarks. A valence clover heavy quark is combined with two im- 
proved staggered light quark propagators to form baryons with different flavors. We are using 
MILC coarse gauge configurations with a lattice spacing of about 0.12 fm. In this preliminary 
investigation, we explore the chiral limit by studying two light sea quark masses, three different 
strange valence quark masses, and nine different light valence quark masses ranging from 0.1 to 
0.4 times the nominal strange quark mass. 
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Table 1: Heavy baryons, flavor contents and operators. J p is total spin parity of the heavy baryon, and s K is 
spin parity of the light quark pair. In the quark content column, / is the light valence quark, s is the strange 
quark, and c is the charm quark. 

1. Introduction 

Heavy baryons have been extensively investigated by experimental and theoretical approaches. 
From experiment, the singly charmed heavy baryon mass spectrum is well known; however, the 
other heavy baryon masses are only crudely known. From lattice QCD, there are several quenched 
calculations for the heavy baryon mass spectrum [1-6], and most results are in fair agreement with 
observed values. In this project, we apply lattice QCD with dynamical sea quarks for the same 
objective. In the future, we hope to study semi-leptonic decays of heavy baryons [0] [Qj. 

We present preliminary results for the mass spectrum of five different singly charmed heavy 
A baryons: A c , E c , S c , E' c , and O c . In this work, we use two different interpolating operators 
studied in Ref. [jl|] and construct two-point functions using the method of Wingate et al. [||] to 
combine staggered propagators for the light valence quarks and a Wilson type (clover) propagator 
for the heavy valence quark. 

2. Construction of two-point functions 

The interpolating operators to describe i + singly heavy baryons are 

where e a bc 1S the Levi-Civita tensor, y/i an d ty% are light valence quark fields for up, down, or 
strange quarks, is the heavy valence quark field, and C is the charge conjugation matrix. As 
can be seen, the spinor indices of the light quark fields are contracted together, so the spinor index 
of the operator comes from ^h- Basically, 0$ is the operator for s n = + , and G„ is for s % = 1 + , 
where s % is spin parity state for light quark pair. Each heavy baryon is obtained by choosing 
different quark flavor combinations, shown in Table |]. 

Since we use staggered fermion propagators for the light quarks, and a Wilson type propagator 
for the heavy quark, we will use the method of Wingate et al. |]] to convert from a staggered 
propagator to a naive quark propagator. They show 

G v (x,y)=£l(x)&(y)G x (x,y) (2.2) 

where, 

&{x) = Y\{Y^ la - (2-3) 
I 1 
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Gy(x,y) is the naive quark propagator, and G x (x,y) is the staggered quark propagator. With this 
basic relationship, we can construct the heavy baryon two-point functions 

C a p(p,t) = £e-^ 5 «(x,0^ 5 /j(0,0)) (2.4) 

X 

= £ e-'Phabce^tr [Gf T (x, 0)Cy 5 Gf (x, 0)Cy 5 \G c ^ (x, 0) . (2.5) 

X 

G\(x,0) and G2{x,0) are naive light quark propagators, Gh{x,0) is the heavy quark propagator, 
and a and /3 are spinor indices. As expected from the property of the interpolating operator, the 
spinor structure of the heavy baryon is determined from the heavy quark. Note that the trace in 
Eq. ^5| is for spinor indices, not for color indices. Now, using Eq. [2.2[ 

Gf (x, 0) = &(x)& (O)Gg' (x, 0) = Gl{x)Gf x (x, 0) (2.6) 

and 

Gf T (x,0) =a T (x)GT x (x,0). (2.7) 

So 

C a p(p,t) = ^e^e^e^tr^ (2.8) 

x 

But the trace can be simplified: 

tr[a J ' {x)Cj 5 a{x)Cy s ] = tr[-l(-lf = -4. (2.9) 

Finally, the two-point function can be expressed as, 

C a p (p, t) = £ e-^e^Sa^A-^GTx (x, 0)G? X (x, 0)G C ^ (x, 0) . (2. 10) 

x 

This two-point function is constructed using the operator 1^5. Similarly, we can derive the two-point 
function with the other operator ff^. 

3. Data analysis 

We use two ensembles of 20 3 x 64 MILC coarse dynamical lattice gauge configurations with 
lattice spacing a f» 0.12 fm. We investigate 385 configurations for the ensemble with ami = 0.01, 
and am s = 0.05, and 458 configurations for the ensemble with ami = 0.02, and am s = 0.05, where 
mi is the light sea quark mass, and m s is the strange sea quark mass. For each configuration, we 
require several propagators for the valence quarks. We compute nine different staggered light quark 
propagators with masses between 0.005 and 0.02, and three staggered strange quark propagators 
with masses 0.024, 0.03, and 0.0415. For the heavy charm quark, we use only one kappa value 
K = 0.122 based on tuning for our heavy-light meson decay constant calculation 

For fitting the baryon propagators, we use correlated least squares fits, and for error estimation, 
we generate 1000 bootstrap samples. Taking into account the periodic boundary condition in time 
for the valence quarks and the staggered phases, the fit model is 

P(t) =Ae-' nt +Ae- m{T -'^ + (-l)'Ar* + (-l)'yfc-*< T -') 

+A Sf e- m * t +A,e- m ^ T ^ + (-1)^^ + (-l)'i«, < T**( I '- f ) ) (3.1) 
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where m is the ground state of positive parity, m has opposite parity, and the values with an asterisk 
denote the corresponding excited states. In practice, most propagators could be fit including contri- 
butions with just two or three particles. In most cases the confidence level of the fit is between 40 
and 60 percent. (Of course, on each ensemble, baryon propagators with different valence masses 
are correlated.) 
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Figure 1: Chiral extrapolation of light valence quark mass on each ensemble. 

To discuss extrapolation to the chiral limit, we present an example for E' c . First, let's look at 
the extrapolation of the light valence quark mass shown in Fig. |]. As we can see, there are three sets 
of data for each ensemble corresponding to each valence strange quark mass. We have performed 
linear extrapolations although there is a very small curvature seen for the lighter sea quark mass 
(Fig. |TJa)). The bursts show the linear extrapolation to the bare light quark mass 0.00148 [1C] 
in lattice units on this ensemble. Next, we interpolate in the strange quark mass (Fig. ^|(a)) and 
extrapolate in the light sea quark mass (Fig. ^(b)). The bare strange quark mass was estimated to 
be 0.039 [10]. In Fig. ^a), we see that linear interpolation in the strange quark mass is perfectly 
adequate. 

After extrapolation in the light sea quark mass, we find the mass of S'. within the partially 
quenched approximation. Alternatively, there is a full QCD point for each ensemble, and we can 
simultaneously extrapolate in light sea and valence masses. Fig. || compares the results of these 
two approaches for four baryons. As we can see, the full QCD results have smaller error bars than 
the partially quenched results. Our multistep partially quenched interpolation and extrapolation 
procedure may be causing larger errors than the simple full QCD extrapolation in light quark mass. 
If we had fit all the data with a partially quenched formula depending on both sea and valence quark 
masses, the partially quenched QCD chiral extrapolation might have given smaller error bars than 
the fit to the limited full QCD data. (We plan to try that in the future.) We present our preliminary 
results here based upon the full QCD extrapolation. 



4. Results 

We present mass splittings in MeV between five different singly heavy baryons in Fig. |](a). 
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Figure 2: Interpolation of strange quark mass (a), and extrapolation of light sea quark mass (b). The bare 
strange quark mass is 0.039, and the bare light sea quark mass is 0.00148. Bursts are the interpolated and 
extrapolated results. 
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Figure 3: Full QCD compared with partially quenched QCD chiral extrapolation. Four different masses of 
singly heavy baryons. Note that Q. c dose not appear, since Q. c does not contain light valence quarks. 



Square data points correspond to the experiment data [11], and crosses with error bars indicate our 
results. The experiment errors are considerably smaller than our own, so we ignore them here. 
Our results agree with the experimental values within errors. In fact, the agreement is so good, we 
wonder if we are overestimating our errors. On the other hand, all the baryon masses are correlated. 

We also present the mass spectrum itself in Fig. |](b). In order to get the physical mass M p i iy 
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Figure 4: Mass splitting (a), and mass spectrum (b) of five different singly charmed heavy baryons with j 
in MeV units. 
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of a heavy baryon, we need to add a constant mass shift A to the calculated mass M ca h i-e. 

M phy =M ca i+A. 



(4.1) 



In this project, we chose the constant mass shift A in a simple way to best match the average of the 
experimental masses. 

A = Average(M exp - M cai ). (4.2) 

We should investigate the momentum dependence of the heavy baryon energy in order to calculate 
the kinetic mass of the state. If we use the dispersion relation, we can take the kinetic mass M^„ as 
the physical mass M p h y [§], 



Afjti, 



\2-[E(p)-E(0)Y 
2[E(p)-E(0)} 



(4.3) 



5. Future study 



We have presented preliminary spectrum results based on dynamical lattice QCD. The study 
can be extended in several ways. We shall investigate singly bottom baryons and doubly charmed 
baryons. Recently, the SELEX collaboration reported the measurement of the doubly charmed 



baryon E+ [12], so it would be timely to compute the doubly charmed baryon spectrum. So far 



we have only studied ^ + states, and we would like to extend the study to | + states. In this way, 
we can investigate the hyperfine splitting problem [||] [§] in heavy baryons. In addition, we can 
increase statistics by using more configurations and ensembles, and additional time sources for our 
propagators. We also need to clarify the taste structure of our baryon operators. 
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